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On the Degree Sequence 
of Random Geometric Digraphs 

Yilun Shang^ 



• Abstract 

' A random geometric digraph G n is constructed by taking {Xi,X 2 , ■ ■ ■ X n } in IR 2 

. independently at random with a common bounded density function. Each vertex X, 

is assigned at random a sector Si of central angle a with inclination Yi, in a circle of 
'. radius r (with vertex Xi as the origin). An arc is present from vertex Xi to Xj, if 

Xj falls in Si. Suppose k is fixed and {k n } is a sequence with 1 <§; /c„ <§; n 1 / 2 , as 
n — ► oo. We prove central limit theorems for fc— and k n — nearest neighbor distance of 
out- and in-degrees in G„. We also show that the degree distribution of this model, 
which varies with the probability distribution of the underlying point processes, can 
be either homogeneous or inhomogeneous. Our work should provide valuable insights 
for alternative mechanisms wrapped in real-world complex networks. 



Keywords: Random geometric graph, Random scaled sector graph, Degree sequence, 
Central limit theorem, De-Poissonization. 



1. Introduction 

■ In random graph theory, degree sequences are among the most elementary and essential 

issues. The random geometric graphs G(X, r) have been well studied in the last decade, see 
the monograph |13j . a short overview [18] and references therein. In order to investigate the 
£f~) • typical vertex degree of G(X n , r n ), Penrose(|14]) defined an empirical process of fc n -nearest 

. neighbor distances in X n , and showed the weak convergence of the finite-dimensional 

distributions of that process, scaled and centered, to a Gaussian limit process. He further 
considered the case k n = k fixed in [13] later. Given a finite point set X S M d and given 
x S X, the ^-nearest neighbor distance means the distance from x to its fc-nearest neighbor 
in X. In the geometric setting, the /c-nearest neighbor distance is often a suitable vehicle 
to deal with degree-related properties of spatial point configurations [12]. It is also closely 
concerned with fc-spacing in statistical testing, which has a number of applications, see 
the book [T7], and is of interest in its own right. 

In this paper we extend the method of Penrose and establish results analogous to the 
ones mentioned above for in-degree and out-degree of random geometric digraphs. Our 
result (Theorem 3) shows that the degree distribution of random geometric digraphs in 
the thermodynamic regime can be either homogeneous or inhomogeneous according to 
different underlying distributions of point processes. In particular, the degree distribution 
is Poisson-like when points are uniformly scattered, reminiscent of that of Erdos-Renyi 
random graphs, see the classic book [1] (Chap. 3); otherwise the degree distribution is 
highly skew (or inhomogeneous), similar with that of many large real- world graphs [7]. 
We also mention that the author was recently able to prove the maximum out /in-degrees 
are almost determined [19] . and this phenomenon has been discovered in Erdos-Renyi 
random graphs pQ. For more discussions, see Section 2.1. 
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Let X n = {X±, X2, • • • , X n }, {Xi} are i. i. d. random variables in R d with distribution 
F having a specified bounded density function /. Let V n = {X\, X2, • • • >^Ar n }, N n ~ 
Poi(n). So V n is a Poisson point process with intensity nf, coupled with X n . Let TL\ be a 
homogeneous Poisson process with intensity A on R d and || • || be I 2 norm on R rf . Standard 
random geometric graphs G(X n ,r n ), G(V n ,r n ) are defined as in [13], that is, G(X n ,r n ) 
(or G{V n ,r n )) has vertex-set X n (or V n ) and an edge X,- L Xj (i ^ j) if | \Xi — Xj\ \ < r n . We 
always assume that r n — > as n — > 00. We now define random geometric digraph models 
to use in this paper as follows: 

Definition 1. (d = 2) Let a E (0, 2n] be fixed. Let y n = {Y\, Y2, • • • , Y n } be i.i.d. random 
variables, taking values in [0,2ti), with density function g. Associate every point Xi E X n 
a sector, which is centered at Xi, with radius r n , amplitude a and elevation Yi with respect 
to the x-axis horizontal direction anticlockwise. This sector is denoted as S(Xi,Yi,r n ). 
We denote by G a (X n ,y n ,r n ) (abbreviated as G n ) the digraph with vertex set X n , and with 
arc (Xi,Xj), i ^ j, present if and only if Xj E S(Xi, Yi, r n ). We can define a Poisson 
version G a (J > n ,yN n -,r n ) (G' n for short) similarly. 

In what follows, we will primarily take g = fbl[o,2Tt)> that is, Yi ~ U[0, 2n). We will 
defer the discussion of the case of d > 3, general probability density function g and even 
other norms to Section 6. Actually, the above model has been first introduced in [4j 
under the name "random scaled sector graph", with d = 2, Euclidean norm and n points 
uniformly distributed in [0, l] 2 . This is an important variant of random geometric graph 
which has been revitalized recently in the context of wireless ad hoc networks, and it is used 
to analyze the performance of wireless sensor networks communicating through optical 
devices or directional antennae, which are significant in mobile communication [TT] . Some 
basic properties and graph-theoretic parameters of this model have also been addressed [H 
EJ[6], using basically combinatorial techniques and discretization. 

The rest of this paper is organized as follows. Section 2 contains the statement of main 
results for d = 2, Yi uniformly distributed. Section 3 discusses the asymptotic results for 
means and degree distribution. In Section 4, we give some moments preparatives for de- 
Poisson. Section 5 includes the proof of main theorems. Section 6 is devoted to higher 
dimension and general probability density function g. 



2. Statement of main results 

We will consider two asymptotic regimes. First, take k n = k E N. Second, let k n — > 00, 
and 

l im _^ = 0. (1) 

rwoo ,/n 



Notice that if we want the sequence {fe n } n >i to converge as n tends to infinity, then the 
above two cases are only choices ( and (HJ is technically needed in the proofs). In the first 
regime, define r n = r n {t) by nr n (t) 2 = t, for t > 0, and in the second, define r n = r n (t) 
by nr n {t) 2 = s(k n + t\/k^), for s > 0, t E R. Here we introduce a tunable parameter t 
to adjust the areas of sectors and t has nothing to do with "time", though we will study 
several random processes with t that evolves. Regulating t allows us to tackle the degree 
sequences in fine details. The reason why we choose such r n is to ensure a non-degenerate 
limit, since nr 2 is a good measure of average degree, see the appendix A of [TO]. We 
emphasize that k n is a crucial parameter which appears in two respects, the scale on 
which the degree distribution tails are studied as well as the scaling for the radius r n . 

Before proceeding, we give some notations to ease statement. For A > 0, let p\{k) := 
P(Poi(X) = k) and for A C Z+ let p x (A) := P{Poi{\) E A). For x E R 2 , let (j), $ be 
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the density and distribution function of standard normal variables. Given i£l 2 , define 
B(x,r) the disk with center x and radius r, and let B n (x,t) := B(x,r n (t)), S n (x,y,t) := 
S(x,y,r n (t)) in both limit regimes. Following Penrose [T3J we set X x := X U {x}, if 
^ is a finite set in IR 2 and x G R 2 . Denote by j^X the number of elements in X and 
X(A) := #(X n A) for A C R 2 . 

We will need some further definitions before we can state our main results. In the rest 
of the paper / max will denote the essential supremum of the probability density function 
/, i.e. / ma x := sup{n : \{x : f{x) > u}\ > 0}. Here and in the rest of the paper | • | 
denotes Lebesgue measure. We assume / max < °° throughout the paper. Next, define the 
level set when k n — > oo as L s := {x G R 2 \sf(x) = §} and let Lf := {x G R 2 \sf(x) > 
^}. We also put a mild restriction on density function /: let R := {x G M 2 |/(a;) > 

0, limsup ^M^Z^fj 3 ^ < with some K < oo, and we always assume F(R) = 1. Let c, c' 

be various positive constants, and the values may change from line to line. 

For Borel set A C M 2 , define ^ ut (t,A), ^° ut {t,A) be the number of vertices in A of 
out-degrees at least k n of G n and G' n respectively. More specifically, 

n 

(,° ut (t,A) = ^l[^(5„(x j ,y i ,t))>fc„+i]n[x i eA] 
i=i 

l [P n (S n (Xi,Yi,t))>k n +l] r\[XiGA] 

i=l 

Similarly, for in-degree we have, 

n 

1 [# {Xj G A'n I Xi G 5„ (X,- , Y 3 - ,t ) } > kn + 1] n [XiGA] 

i=l 

N„ 

£™{t,A) = 'Y^\#{XjeT n \x i eS n {x j ,Y j ,t)}>k n +i][\[x i eA] 

8=1 

Notice for the case k n — > oo, s is suppressed in the above expressions. Also, let : = 

£° ui (i,IR 2 ) etc. for convenience. 

The following two lemmas are intermediate steps to prove Theorem 1 and 2. We 
choose to state them without proof due to the limitation of space and they can be treated 
in parallel with Theorem 4.12 and 4.13 in [13] through a dependency graph argument. 
Lemma 1. Suppose that k n = k is fixed, and that A is a Borel set in R 2 . The finite- 
dimensional distributions of the process 

n- 1 2[C° Ut (t,A)-EC out (t,A)} , t>0 

converge to those of a centered Gaussian process (^^(t, A),t > 0) with covariance 
E[^\t,A)^ u \u,A)\ given by 

J P%tf(x)([k,°o))f(x)dx 

+ / / / / ^(^/(^i),yi,y 2 )/ 2 (xi)dzdx 1 dy 1 dy 2 
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with 

^ t (z,X,y 1 ,y 2 ) =P({H z x (S(0,y 1 ,t 1 *)) > k} n {H x (S(z,y 2 , U ^) > k}) 
-P(Hx(S(0,yi,V)) > k)P(H x (S(z,y 2 ,J)) > k) 

The finite- dimensional distributions of the process 

n-H^(t,A)-EC n (t,A)] , t>0. 

converge to those of a centered Gaussian process ((™(t,A),t > 0) with covariance 
E[Z'™(t,A)^u,A)} given by 



with 



p« tf{x) ([k,oo))f(x)dx + I l_y^(z,—f(x 1 ))f%x 1 )dzdx 1 



#S(*,a) =P(M(B(o,il)) > k}n{n° x (B(z,J)) > k}) 

-P(W A (fl(0,t*)) > k)P(H x (B{z,J)) > k) 

Let W denote homogeneous white noise of intensity tt^ 1 on M?, that is, a centered Gaus- 
sian process indexed by bounded Borel sets in M 2 , with covariance Cov(W(A), W{B)) = 
-|j4n.B|, where | • | as mentioned before is Lebesgue measure. Also, let W' denote homo- 
geneous white noise of intensity — . 

Lemma 2. Suppose that k n — ► oo ; that |7P holds, and that A is a Borel set in M. 2 . Let 
s > and suppose F(A n L s ) > 0. The finite- dimensional distributions of the process 
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(nk n )-nC° U \t,A)-EC out (t,A)] , iel 

converge to those of a centered Gaussian process (Z£ ut (t,A),t e R) with covariance 
E[^ ut (t,A)^ ut (u,A)} given by 

\L S n A\ f 2n f 2n f 

J„ n \2 / / / o Cov ( 1 [W'(S(0,y 1 ,l))<t]' 1 [W(5( Z ,y 2 ,l))< U ])d2d 2 /idy 2 

s(na) J J j R 2 
The finite- dimensional distributions of the process 

(nk n rH^n n (t,A)-E^(t,A)} , i£R 

converge to those of a centered Gaussian process (£™(t,A),t £ R) with covariance 
E[^(t,A)^(u,A)) given by 



4- |L s nA| 

sa 2 



Cov(l[ W ( S (o,i))< t ] , l[w(B(z,i))<u])dz. 



Now we are ready to state our main results. 

Theorem 1. Suppose that k n = k is fixed. The finite- dimensional distributions of the 
process 

n^[C ut (t)-EC u \t)] , t>0 
converge to those of a centered Gaussian process {S^it)^ > 0) with 

m^m^iu)] = E[te*{t)&*(u)] - h(t)h(u), 
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where 



h(t) = J {p S t f( x )(k - + p ¥ f( x) ([k, oo))} /(x)ds (2) 

T/ie above result also holds in the case where the superscripts 'out' are replaced by 'in' 
everywhere. 

Theorem 2. Suppose that k n — > oo, and (QP holds. Let s > and suppose F{L S ) > 0. 
The finite- dimensional distributions of the process 



i 



(nk n )-HC ut (t)-EC ut (t)] , tern 

converge to those of a centered Gaussian process £ R) wif/i 

MC'WCV)] = - g(t)g(u), 

where g(t) = 4>(t)F(L s ). 

The above result also holds in the case where the superscripts 'out' are replaced by 'in' 
everywhere. 

To deal with the degree distribution, let 77°"* (i, A) and rf™(t,A) be the number of 
vertices in A of out-degree and in-degree k fixed in G n respectively. 

Theorem 3. Suppose A is a Borel set in M? and a > n. If either k n = k fixed, or k n — > 00 
and n~ A^lnn —* 0, then 

lim n- x C*(*, A) - ^[« _1 C*(*> 4)] = a.e. (3) 

Moreover, 



lim n X (M) = / P%tf{x)(k)f{x)dx a.e. (4) 

The above result also holds in the case where the superscripts 'out' are replaced by 'in' 
everywhere. 

2.1 Discussion of Theorem 3. 

We take expectation on both sides of (J1J), and let p(k) := £ , lim n ^ 00 n _1 ?7° u *(t, IR 2 ), so 
the out-/in-degree distribution of G a (X n , y n , r n (t)), where nr n {t) 2 = t, is 

p(k) = ^- [ e-^ tfix) f(x) k+1 dx, A;GNU{0} (5) 

k\ JR2 

If we take the uniform density function f{x) = l^ 01 p(x) in (|5|), then we see that 
p(k) = e~ 2"*(^/j) fc /A:!, k > 0; that is, the degree distribution is Poi(^t). 

If we take the standard multivariate normal density function f(x) := f(x±,X2) = 
(l/2n)e~( Xl+X2 ^ 2 , then through the polar coordinate transformation and integration by 
parts, we obtain p(k) = (An/at) — e - "*/ 4 " ^( c =0 (Qi/4Tx) 4_1 /i!, k > 0. It is easy to see 
that p{k) — > as k — > 00; and furthermore, since p(0) = (4n/ai)(l - e~ at / 4n ), p(0) -> 1 
as i — > and p(0) — > as i — > 00. These observations allow us presumably adjust the 
parameter i to get different skew degree distributions especially for small k. However, the 
degree distribution in ([5|) has a light tail in contrast to the power law distributions [7] 
because of the fast decay as k tends to infinity. To be precise, by © and Stirling formula, 

( ft/max) fc f ( , _, , „ (ate/ max ) fc p 
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for any (3 > as k — > oo. 

On the other hand, if we want to find a suitable density function / for a given prob- 
ability distribution p(k) satisfying p{k) > and YlT=oP^) = 1> then we simply solve 
the equation ([5]), which is the first kind nonlinear singular Fredholm integral equation [3]. 
However, only approximation solutions of this kind of equations may be obtained by using 
iterative methods and the existence of solution is not known in general. 



3. Proof of means and degree distribution 

Proposition 1. (out-degree) Suppose A C M? is a Borel set. If k n = k is fixed, then 



limn~ 1 E[C ut (t,A)]= / p« tf(x) ([k,oo))f(x)dx 



If k n — > oo, and (Op holds, th- 



en 



lim n~ l E[i° n ut {t, A)] = F(L+ ni) + ${t)F{L s n A) 

n— too 



(6) 



(7) 



Proof. Let p n (x,y,t) = F(S n (x,y,t)). Then 



E[t™%A)) = ^j* U J a P[B 



in(n- l,p n (x,y,t)) > k n ]f(x)dxdy 



(8) 



Suppose k n is fixed, and x £ R, then / is continuous at x and np n (x,y,t) — ► ^tf(x) 
by mean-value theorem of integrals. Therefore P[Bin{n — l,p n (x,y,t)) > k] tends to 
p^tf(x)([k, oo))- Then ([6]) holds by (|8]) and dominated convergence theorem. 

Suppose k n — > oo, ([I]) holds and x £ R, then np n {x,y,t) ~ n^r^f(x) ~ s~f(x)k n , and 
by Chernoff bounds (see e.g. [5]), P[Bin(n — l,p n (x,y,t)) > k n ] tends to 1, if s/(x) > ^; 
and tends to 0, if sf(x) < ^. Then for i£fin L s , 



np n {x,y,t) = n-rlf(x) + n 



S n [x,y,t) 



(f(z) - f(x))dz 



= k n + tkl + 6(n(A: n /n) 3 / 2 ) 

Hence, by ([1]), 

i i 

np n (x,y,t) = k n + tk£ + o(k£), xeRDL s 
Then let p n = p n (x,y,t), by DeMoivre-Laplace limit theorem and Q, we have 

P[Bin(n - l,p n ) > k n ] 

= P 



0) 



Bin(n - l,p n ) - EBinj n - l,p n ) k n - (n - l)p ri 



y/np^ 

So ([7]) follows from ([8]) by dominated convergence theorem. □ 

Proposition 2. (in-degree) The same results hold when replace superscripts "out" by "in" 
in Proposition 1. 

Proof. Let q n (x,t) = ^ ■ F(B n (x,t)). Then 

E[C(t,A)} = n [ P[Bin(n - l,q n (x,t)) > k n ]f(x)dx. 
J A 
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From Palm theory, similarly we have 



E[& n (t,A)} =n [ P[Poi(nq n (x,t)) > k n ]f(x)dx. 
J A 

The remain proof is in a similar spirit with that of Proposition 1. Hence we omit it. □ 
We remark here that Proposition 1 and 2 still hold for corresponding Poisson case. 

Proof of Theorem 3. Define a a filtration: Tq = {0,(7}, and for 1 < i < n, Ti = 

a{(X 1 ,Y 1 ),(X 2 ,Y 2 ),--- ,(X h Yi)}. 

For out-degree, £*(t, A ) ~ A)] = £? =1 M%* , with Mf£ = E[^(t, A)\Fj\ - 

E[£° ut (t, A)\Fi-i]. Let Z$(t, A) be the number of vertices in A of G(X n+1 \{Xi},y n+1 \{Y} 
,r n ) having out-degree at least k n . Thereby, M™' = E[^ ut (t,A) - $$(t,A)\Fi]. 

We now claim that: For finite set <YCM 2 and x G X, there are at most 8k points z £ X 
having x as their (< k) — th nearest neighbor, for any k £ N. Here x is the k — th nearest 
neighbor of z in X means if we order quantities \\\w — z\\ : w £ X\{z}} increasingly, then 
\\x — z\\ will be the k — th item in this sequence. Proof. We take a cone with vertex x, 
central angle ti/4. It's easy to see that there are at most k n points of X having x as their 
(< k) — th nearest neighbor, since we may look for these points from near to far. The 
claim follows since the plane is covered by 8 such cones. 

Therefore, 

\C u \t, A) - £?(t, A) | < |C*(*, A) - A) | + Igftfe A) - $g(t, A)\ 

< (8k n + 1) + (8k n + 1) < 18k n , 

where let £n+i(*> A) denote the number of vertices in A of out-degrees at least k n of 
G(X n+ x,y n+ x,r n ). Then |M?jf| < 18k n . For e > 0, by Azuma inequality, see e.g. [2], 

P[|C*(M) - E[g*(t,A)]\ > en] < 2e - £2 « 2 /648n/4_ 

By Borel-Cantelli Lemma, ([3]) follows. The in-degree case can be proved similarly. 
To prove (jlj), we notice 

n n 

r]° ut (t, A) = } j \x n (s n (Xi,Yi,t))>k n +i] CiiXieA] ~ X] 1 [A' n (5„(x Ii y,,t))>fc„+2] Dlx^A] 

i=l i=l 

and by (|3|) and the proof of Proposition 1, the result follows immediately. The in-degree 
case also follows similarly. □ 



4. Some moments for de-Poissonization 

In this section we will develop some moments for non-Poisson case in the limit regime 
k n — > do, which is crucial to de- Poisson Lemma 1 and 2. 
For n, m £ N, set 

m 

r m"n(*) := l [X m (S n (X l ,Y l ,t)\{X 1 })>k n ] 
i=l 

and 

m 

Tm^nit) '■= ^2 1 [#{X- i eX m \{X 1 }\X 1 &S n {X :j ,Y 3 ,t)}>k n ] 
i=l 
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Then we see T^(t) = ^(t),T^ >n (t) = and I*» ft (t) = &(t),T™ >n (t) = Q n (t). 

Set D% n (t) := T™t 1>n (t) - I™* (t), then D^ n (t) = D^ n (t) + D%{t), where 

rn 

D m%(f) = 1 [X m (S n (X l ,Y l ,t)\{X l })=k n ~l]n{X m+1 £S n (X l ,Y>,t)] 
i=l 

Dmfnit) = l [X m (S n (X m+1 ,Y m+1 ,t))>k n ] 
Set D™ n (t) := T- +ljn (t) - I™ n (t), then ££ n (i) = £# >n (t) + D™ n (t), where 



m,n 

(*) = £ 



j ,*) }=fcn - 1] n [XiGS„ (X m+1 ,Ym+l ,*)] 



i=i 



D^nit) - 1 [#{x J eA' m |x m+1 e5„(x j ,y J ,t)}>fc, l ] 

We denote binomial probability (3 np (k) := P(Bin(n,p) = k). The next lemma will be 
repeatedly used in this section, see [TBI [T4"]. 

Lemma 3. (a,) Suppose n, k £ N loit/i k < n. Then f3 n ,p(k) is maximized over p £ (0, 1) &?/ 
setting p = k/n, and pf3 n , p (k) is maximized over p £ (0,1) fry setting p = (k + l)/(n + 1). 

(b) Suppose {j n }n>i is a sequence of integers satisfying j n — ► oo and (j n /n) — > as 
n — ► oo. Suppose t 6 R and {p n }n>i «s a sequence in (0, 1) satisfying {jn—np n )/(np n ) l l 2 — > 
t as n — > oo. T/ten 

Jn 2 Pn,p n (jn) ~> CIS 71 > OO. 

Lemma 4. Suppose k n — > oo and (QP holds. Then 

lim sup |^^»t w _^ )i?(Ls)|=0 

{m||m— n|<ra 2 / 3 } 

TTie same formula holds when replace superscript "out" by "in". 

Proof. Take {m n } n >i with |m n — n| < n 2 / 3 . 
For out-degree, we have 

k-^EDZUt) = T~ 1^ I m n k~ 1 / 2 P(X mn ^ 1 (S n (x,y,t)) = k n -l) 
ln Jo Jr 2 

■F(S n (x,y,t))F(dx)dy. (10) 

Let x G RC\L S , then A' mn _i(5 n (x, y, t)) is binomial with parameters m n —l and F(S n (x, y, t)), 
and by Q, (P) the mean is 

m n F(S n (x,y,t)) = (l + 0( n - l ^))(k n + tk^ 2 + o(k l J 2 )) 

= k n + tk\[ 2 + o^ 2 ), x £ Rn L s (11) 

By Lemma 3, 

lim k*/ 2 P(X mn - 1 (S n (x,y,t)) = k n -l)) = <l>(t), x £ R Pi L s 

n—>oo 

Also, by Chernoff bounds and Proposition 1, 

lim $/ 2 P(?C mn - 1 {S n {x,y,t)) = k n -l)) = 0, x £ R\L S 
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Hence for x G R, the integrand on the right hand side of (fTUj) tends to 

1/2 

Also, by Lemma 3, (m n /k n )F(S n (x, y, t)) and k n sup 0<p<1 0m n -i,p\k n — 1) are uniformly- 
bounded. So, fc n ED™* n (t) tends to <p(t)F(L s ) by dominated convergence theorem. 

Since < D™^ n (t) < 1, kn 1/2 ED^ n (t) tends to 0. The first part of the lemma then 
follows. 

For in-degree, we first introduce some notations. Let / := and for Borel set 

A C M 2 , let ^(.A) ~ Biti(u,F(i)), where F(A) := J A f{x)dx. 
Consequently, we have 

k- 1 / 2 ED% n>n (t)= [ m n k~ 1 / 2 P(X m ^ 1 (B n (x,t)) = k n -l)F(B n (x,t))F(dx). 

Let x G R(~)L S , as mentioned above, X mn -i(B n (x, t)) is binomial with parameters m n — 1 
and F(B n (x,t)), and by Proposition 2 and ([!]) the mean is 



n F(B n (x, t)) = k n + tk\l 2 + o{k]j 2 ), x £ R H L s 



m 

By using Lemma 3 and Proposition 2, we can conclude the proof in a similar manner 
with the out-degree case. □ 

Lemma 5. Suppose k n — > oo and ([7J) holds. Then 

lim sup \tfEBff(t)D°g n (u) - m<P(u)F(L s ) 2 \ = 

n—n 2 / 3 <l<m<n+n' 2 / 3 

The same formula holds when replace superscripts "out" by "in". 

Proof. Let I < m. 

For out degree, we have 

l m 

EDff(t)D% n (u) = ^^P[{AKS„(Xi,^)) = K} n {X m {S n (X v Y v u)) = k n } 

i=l j=l 

n{x l+1 e s n {Xi,Y h t)} n {x m+l g s^x^u)}] 

r-2n r 2n 

9n,l,m 

(xi,yi,X2,y2)F(dx 1 )F(dx2)dyidy 2 



4ti 2 

l{m-l) 



Jo 



2n r 2n 

11 i 

An 2 



9n,Lm( x i^yi^ x 2,y2)F(dx 1 )F(dx 2 )dy 1 dy 2 



I 

'2n 



2n 







Jo 

9n,Lmi x ^yi) F ( dx l) d yi ( 12 ) 



where, 

9n,l,m{ x l,yi, x 2,y2) ■= P\{^-2 ( S n( x l , 2/1 , *)) = K - 1} f] {AVi G S^l, 2/1, *)} 

^{ X m-2( S n( x 2,y2,u)) = k n - 1} H {X m _i G 5 n (x 2 , 2/2, «)}] 

9n,i, m ( x i,yi, x 2,y2) ■= P[{Xl-i (S n ( x l , 2/1, *)) = A; n - 1} n {X/ G S^l^l,*)} 

n{^ 1 _ 2 (5' n (x 2 ,y 2 ,'u)) = fc n - 1} H {X m _i G 5 n (x 2 , J/2, «)}] 

S^mO^iS/l) := ^[{^-1(^(^1,2/1,*)) = /c„ - 1} n {X; G S n (xi,yx,t)} 

r\{Xm-i(S n (xi,yi,u)) = k n - 1} n {Jf m G SW0ei,2/i,«)}] 
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Take xi,X2 G R, x\ ^ X2 and 2/1,2/2 G [0,2n). Take {Z n }n>i and {m n } n >i with 
n — n 2 / 3 < l n < m n < n + n 2 / 3 . Then asn-> 00, 

r P(^-i e S n (x 1>yi ,t)) - — /(xi), (13) 

T7 CA/ 

— P(X mn _x G S n (x 2 , y 2 , u)) -> — /(x 2 ). (14) 

Since 



^(^,-20^(^2, 2/2, «)) = fc n - l|Xj n _l G S n (xi,yx,t)) 

P(Xm n -2(S n (x 2 ,y2,u)) = k n ) • P{X ln _i G S" n (xi,yi,t)) 



))(^n), 



P(X in _! G S n (x uyi ,t)) 

by Lemma 3 and (jlip . we obtain 

fcy 2 P(^_ 2 (5 n (x 2 ,y 2 ,u)) = 1|X, B _! G SnCxi, !/!,*)) - 0(u), x 2 G L s (15) 

Let Xi,x 2 G M 2 with X2 G" B(x\,r n {t) + r n (u)), so B n (x\,t) C\ B n (x2,u) = 0. If 
^ mn _ 2 (5' n (x2,y2,M)) = fc n - 1, then A/ n _ 2 (S' n (x2, 2/2, it)) = j, for some < j < k n - 
1. Given <^ n _2(<S n (x2, 2/2) it)) = j, the conditional distribution of ^i n -2(S n (xi, yi, t)) is 
binomial with parameter l n — 2 — j and F(S n (xi,yi,t))/(l — F(S n (x2,y 2 ,u))). For all such 
j, if also xi G L s then by (fTTj) . (pQ) the mean of this distribution is 

(i„-2-i)F(S„(z,,, l ( )) _ + t4;/ 2 + 0(t a/2 ))(1 + 0( ^ )} 
l-F(iS n (x 2 ,y2,«)) n 

Therefore for x\ G L s and £2 7^ xi, by Lemma 3 we have 

feV^t^CS'n^.yi,*)) = fcn - l|{^C_ 2 (5n(x 2 ,|tt,«)) = ^ - 1} 

n{Xi„_! G 5„(xi,2/i,i)}] ->^(t). 



Combining this with (|T3j) . (fill) and (jl5jl . we get 

(n 2 /k n )g n j n!mn (x 1 ,y 1 ,x 2 ,y2) -> 4>(t)cf)(u), x 1 ,x 2 £ Rn L s ,x 1 ^ x 2 (16) 
On the other hand, by Chernoff bounds, 

k nP[X^ 2 {S n (x2,y2,u)) = k n - l|X Jn _i G 5„(a;i,yi,t)] 0, x 2 G -R\L S 

and 

^[^^(^(xi, yi, i)) = fc n - 1] -»■ 0, X! G ii\L s 
Combing these with (|13|) and (|14p . we have 

{n 2 /k n )g nt i ntmn (xi,y 1 ,x 2 ,y2) 0, (x x , x 2 ) G (i? x R)\(L S x L s ) (17) 

If x 2 P(xi,r n (t)+r n (n)), setting pi = F(5 n (x 2 , 2/2, it)) andp 2 = F(5 n (zi,yi,t))/(1- 
pi), we get 

£n,/„,m„ (^1^1,^2,2/2) < max PlP2Ann-2.p1 (&n - l)Pl n -2-j,p 2 (k n - 1) 

0<j<fc n -l 
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Whence by Lemma 3 and Stirling formula, there exists a constant c such that 

-I k 

9n,l n ,m n {xi,yi,X2,y2) < c{k n 2 • — ) 2 = ck n /n 2 . 

n 

Then by (|16p . (|17|) and dominated convergence theorem, we get 

Um L re / | TT 2 l l l I 9n,l n ,m n ( x l,yi^2,y2) 
n^oo fc n -4Tt J J Q J M 2 JM?\B(xi,r n (t)+r n (u)) 

F(dx2)F(dx 1 )dy 1 dy 2 ] = <t>(t)<t>{u)F(L s ) 2 . (18) 



Also, by C3P and 

9n,l n ,m n ( x l^yi^ x 2,y2) < P { X l*-2 ( S n ( x l , 2/1 , t) ) = fc n - 1)F(5„ (x X , J/l , t) )F(S n (x 2 , 2/2 , It)) 



n 

Since F(B(xi,r n (t) + r n (u))) < c • (k n /n) for some constant c, by (H|), 

,2 /•2-n r 2n 



( , 2 ) I I I I 9n,l n ,m n ( x l,yi,X2,y2) 

tin ' 4ti Jo Jo j r2 JB(a; lirn ( f ) +r „(„)) 



F(dx 2 )F(dx 1 )d 2/1 d 2/2 < t /(_)(-!L)(-rL_) - 



3/2 
n 



Thus (fTBj) holds with the region of integration modified to [0, 2tc) x [0, 2tc) xl 2 xR 2 . The 
asymptotic results for g' n h m are just the same. Also, by similar arguments there is a 
constant c such that 

l n ■ sup g" ln (xi,yi) < cnk-^ikn/n) 2 -> 0. 

Vl€[0,2n) 

Hence flT2|) yields 

k-'ED^D^u) -> <j>{t)<l>{u)F(L 8 ) 2 . 

What remains to show is that the above formula still holds when is replaced by 

D£f n ; is replaced by D££ n . We have < Dff n (t) < 1, < < 1. 

By the proof of Lemma 4, EDf£ n (t) = 0(k n /2 ) and ED%f n n (u) = 0(k n /2 ). Therefore 

E[D^ n (t)D^ n (u)}, E[b°^{t)D™l n {u)] and are all 0{k n ' 2 ). The 

first part of this lemma whereby follows. 

The proof for in-degree case parallels to the above approach and we leave it as an 
exercise for the reader. □ 

Lemma 6. Suppose k n — > oo and (T7]) holds. Let t, u G R. Then 

limsup (k~ 3 / 2 ■ sup E[D^ n (t) 2 }) < oo. 

n *°° {m||m— n|<n 2 / 3 } 

The same formula holds when replace superscript "out" by "in". 
Proof. Take {m n }n>i satisfying \m n — n\ < n 2 / 3 . 
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For out-degree, by (|12|) with I = m = m n , t = u, 

E[Dmln{t?\ = mn{m A n ~ 1} [ [ 5„, m „, m „(x 1 ,y 1 ,x 2 ,y 2 )F(dx 1 )F(dx 2 )dy 1 dy 2 

4ti^ J Jo J r2 J r2 

By Lemma 3, there is a constant c such that 

gn,mn,m n ( x U Vl,X2, V2) < sup [pmax{/3 m „, p (fc„ - I), (3 mntP (k n - 2)}] < ( — )k~ 1/2 . 

o<p<i n 

Also, g n , mn ,m n (xi,yi,x 2 ,y2) = unless x 2 G 5(xi, 2r n (t)). Whence 

/ 9n,m n ,mn 

(xi,yi,x 2 ,y 2 )F(dx±)F(dx 2 )dyidy 2 < — 5— 



10 Jo 



By Lemma 4, ED^ n (t) = 0(k l J 2 ) and m n (m n -l) = 0(n 2 ). So E[D^ >fl (t) 2 ] = 0(kl /2 ). 
The first part of this lemma then follows, by noting < D^n n(^) — 1- 

For in-degree, the same argument may be applied. Thus we conclude the proof. □ 



5. Proof of central limit theorems 

To prove Theorem 1 and 2, we will employ useful de-Poisson techniques given in [9], 
[15] and later generalized in [131 US]- We will also need Cramer- Wold device, see e.g. |5]. 
Now we are in position to prove our main results. 

Proof of Theorem 1. Let M G N, B = {h,--- ,b M ) G R M , T = (h,--- ,t M ) G (0,oo) M . 
For out-degree, ^ctfjC [0,2ti) with card(Ar) = card(3 ; ), set 

M 

H (X,y):=Y^ Yl 6il [*(S(x lW ,^ 3 ))>fc n +l] 
*=i (a !) v)e(Af,y) 

and let H n (X,y) = H (n^ 2 X,y). (x,y) G K 2 x [0,2ti) C M 3 . Set £° ut {T,B,A) := 
E^^^A) and Var(C;° u '(T,S,A)) := a' out (T, B, A), we have H n (V n ,y Nn ) = 
£° ut {T,B,R 2 ), and what's more, (V n ,yN n ) is a 3— dimensional Poisson process, which 
may be coupled with (X n ,y n ) in the same way as V n does with X n . By Lemma 1, 

n- l l 2 {H n {V n ,y Nn )-EH n {V n , y Nn )) AA(0, a' out (T, B, M 2 )). Let W A be a 3-dimensional 
homogeneous Poisson process and denote point (x,y) G R 2 x M, W A := (H^\H { ^) 

with x G H A 1} ,y G W A 2) . Next, we say floO^jD^) is strongly stabilizing on Ti.\ if there 
are a.s. finite random variables T and A(W A ) such that with probability 1, A (A) = 
A(H\) for all finite A := (A 1 ,A 2 ) C R 2 x [0, 2tc) with card(Ai) = card(A 2 ), satisfying 
An (5(0, T) x [0,2ti)) = H x n (5(0, T) x [0,2tx)). Here, A(W A ) := H (H° X ) - H (H X ). 



12 



Thus, Hq is strongly stable since it has finite range. We have 
E[A(H X )) = E[H (H° X ) - H (H X )} 

M 

= E [Y.H Yl 1 [H{ 1)fi (S(x,y,ty 2 ))>k+l] +1 {H{ 1) - a (Smt] /2 ))>k+l]} 

t=l (x, y )en x 

M 

E \H£Hs(x,v,t\' a ))>k+l])] 
i=l (x,y)eH x 

M 

= EY J H\ Poi( 2n\-^)>k\+ S 1 [W«(5(x l2/ ,^'' 2 ))=fc]) 
i=l ^,v)en x 

M 

= J2bi(p Xn0(ti {[k,oo))+\2n-^(k-l)). 

i=i 

By © and the Cox process H^ X ,Y) with Y) := (1/2ti)/(X), we have £[A(7^( X ,y))] 
= X^i Set *max = max{ti, • • • ,t M }, we have \H n (X m ,y m )\ < mj^fii \k\ a nd 

M 

\H n (X m+ i,y m+ i) — H n (X m ,y m )\ < (/' J bi) ■ [#{Xi £ X m \X m+ i G S n (Xi,Yi,t ma _ x )} + l] 

i=l 

which is stochastically dominated by c • [Bin(ffl, /ma X Tir n (t m ax) 2 ) + 1] having a uniformly 
bounded fourth moment when m < 2n. Therefore by a simple variant of Theorem 
2.16([13j) to a marked point process [16] (in particular the translation-invariance of Hq(X, y) 
is only required for X), n- l l 2 (H n (X n ,y n ) - EH n (X n ,y n )) M(0,T* ut ) with r 2 out : = 
cj' ou \T, B, R 2 ) - (E[A(H^x,Y))]) 2 - The first P art of the theorem then follows by Cramer- 
Wold device. 

For in-degree, let X C M 2 ,3^ C [0, 2tx) with card(A') = card(^), and the elements 
(x,y) G (X,y) be ordered pairs. Reset 

M 

H {X,y) :=J2J2 bil [#{^xW£S(x,y^ 2 )}>k n +i] 
i=i x'ex 

and let H n (X,y) = H^n^X, y). Set £ n (T,B,A) := Em=l b m & n (t m ,A) and 
Var(#"(T,B, A)) := a' in (T, B, A), we have H n (V n ,y Nn ) = & n (T,B,R 2 ), and (V n ,y N J 
is a 3— dimensional Poisson process coupled with (X n ,y n ). By Lemma 1, n~ l l 2 {H n {Vri, JV n ) 
- EH n (V n ,y Nn )) 7^(0,cr' in (T,5,]R 2 )). Also, H is strongly stable. Let H x be a 



13 



3— dimensional homogeneous Poisson process and TL\ := (T~t\ ,7~C\ ) as above. Then 
E[A(Hx)] = E[H (H° x )-H (Hx)} 



M 



- E[^2 bi( 1 [#{^n{ lh0 \x'es(x,y,tl /2 )}>k+i] + 1 l#{^H[ 1) ' \oeS(x,y,t] /2 )}>k+i]) 



- 1 x'eH^ 



M 



^Jbj( 1 [#{a-GW^ 1) |x'g5(a:,s/,t 1 1/2 )}>fc+l])] 



=1 



- E ^2H 1 [Poi(2nX-^)>k] + X] 1 [#{xGWl 1) |a''e5'( a; ,j / 4 /2 )}=fc]) 

i=1 K / ew[ 1) n5(o,o,tV s ) 



i=l 

The remain proof is similar with the out-degree case. □ 

Proof of Theorem 2. Let T and B £ R M . 

For out-degree, <YcR 2 ,K [0, 2tc) with card(^) = card(3 7 ), set 

M 

i=l (a;,H)e(Ar,y) 

By Lemma 2, we have H n (P n ,y N J = kn l/2 ^ ut {T : B,M?) and n- x l 2 {H n {V n , 3VJ - 
EH n (V n ,y N J) AA(0,<r w (T,£,]R 2 )). Set a := Y^=xbi<f>{U)F{L s ), and fl£f n := 
" ^n(^ m ,^ m ). Then = k- ll2 Y!tihD°^ n {ti). By Lemma 4, 5 

and 6, we have 

lim ( sup \EKti -a\) =0 



n— n 2 / 3 <m<n+n 2 / 3 



lim ( sup |£[/2^ t n i^J-a 2 |)=0 

n— n z ' J <m<m' <n+n' 1 / ° 

and 

lim (n- 1 / 2 sup E[(R° n %) 2 }) = 0. 

respectively. Also \H n (X m ,y m )\ < mj]^ Then Theorem2 . 1 2 ( [IB] ) implies n -1 / 2 

• (ffn^,^) — EH n (X n ,y n )) M(0,a out (T, B)), with cr mt (T,B) := <r w (T, 5, M 2 ) - 
a 2 . Hence, a OM '(T, B) = Var h&ifi)- The first part of the theorem then follows by 
Cramer- Wold device. 
For in-degree, let 

M 

H n {X,y) := k~ 1/2 ^2 X h i l [#{^X\x'£S n (x,y,U)}>k n +l] 

j=l x'ex 

We then argue likewise to complete the proof. □ 



6. Further discussion and remarks 
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In the above sections, we consider d = 2 and Yi uniformly distributed. A natural 
generalization is to consider higher dimensions. For example, for d = 3, instead of a 
sector with amplitude a, we have to consider a spherical sector SS(X, Y, Z, r) which is 
the region bounded by a cone with vertex X, central angle a and a sphere with center X 
and radius r. We take X as the origin and build the standard right-handed coordinate 
system. Let the chief axis of the cone be a ray I, project I onto xOy— plane, and call it V . 
Let Y be the angle between positive x— axis and V and Z + (a/2) be the angle between 
I and V . Y, Z £ [0, 2tx). Consequently, the formal definition of this "random spherical 
sector graph" is easily stated. If Y and Z have uniform distribution, and instead of 

2 

condition ([I]) we assume k n /n d + 2 tends to and modify the definitions of r n (t) accordingly, 
then analogous results corresponding to those appeared in above sections can be derived. 
Actually, we have for example, £,° ut (t,A) = 1 [A' n (S5„(x l ,y l ,z l ,t))>fc n +i] HPQeA] • Let 

p n (x, y, z, t) = F(SS n (x, y, z, t)), then 

,-2n ,-2n r 

£[C'(M)] = ^ / J Q J A P[Bin(n-l,p n (x,y,z,t))>k n ]f(x)dxdydz. 
Also, Cn(t,A) = Y£=i 1 [#{x J eA' n |x ! e5S„(x J ,y j ,z j ,i)}>fc Il +i]n[^eA]- Let q n (x,t) = F(B n (x,t)) 

.( l-°°|("/2) ) >then 

E[C(t,A)} = n I P[Bin(n - l,q n (x,t)) > k n ]f(x)dx. 
J A 

Another direction to investigate is to consider probability density g of Y other than 
the uniform density. Suppose EY < oo. For out-degree case, we may proceed smoothly by 
similar argument, whereas for in-degree the story is different. Say, we consider in-degree of 
a vertex u. Suppose \\u — v\\ < r. Since the inclination of sector S v now is not uniformly at 
random (as we now consider a general density g), we will have distinct thinning probability 
for different v. Moreover, the probability of vertex u lying in the sector S v essentially relies 
on not only the distance between them but also the position of both vertices u and v. Then 
the computation is inevitably involved and the above de-Poisson technique is no longer 
valid. 

We mention that the model is less interesting when using other non-Euclidean norm in 
application viewpoint. It is easy to see when d = 2, if we take l p (1 < p < oo) norm, and 
a = tx/2, Tt, 3tc/2 or 2tc, the above results still hold, due to the symmetry of the coordinate 
vectors under such norm. 
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